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We investigate hadron production in deep inelastic lepton-nucleus scattering in the 
kinematic regime of the HERMES experiment. Our calculations are carried out in 
the framework of a BUU transport model which contains the Lund event generators 
PYTHIA and FRITIOF for the simulation of high-energy elementary interactions. 
For the first time we consistently use the complete four-dimensional information of 
the Lund string break up vertices as input for our transport theoretical studies of 
(pre)hadronic final state interactions. We compare our results with experimental 
HERMES data on charged hadron attenuation. 



1. Introduction 

Hadron production in deep inelastic lepton-nucleus scattering (DIS) pro- 
vides an ideal tool to investigate the space-time evolution of hadron 
formation 1 . The nuclear target can be viewed as a kind of 'micro-detector' 
that is located directly behind the virtual photon-nucleon interaction ver- 
tex. It allows us to study the interactions of the reaction products with 
the surrounding nuclear environment on a length scale that is set by the 
size of the target nucleus. By comparison with hadron production on a 
deuterium target one can draw conclusions about the space-time picture 
of hadronization. The latter information is crucial for the interpretation of 
jet quenching in ultra-relativistic heavy ion-collisions at RHIC as a possible 
signature for the creation of a deconfined quark- gluon plasma phase 2 . 

In the recent past the HERMES collaboration at DESY has started 
an extensive experimental investigation of hadron attenuation in DIS on 
various gas targets 3 and a similar experiment is currently performed at 
Jefferson Lab 4 . Furthermore, hadron attenuation in nuclear DIS will be 
subject to investigation after the 12 GeV upgrade at Jefferson Lab and at 
a possible future electron-ion collider 5 . In Refs. 6 the observed attenuation 
of high energy hadrons at HERMES has been interpreted as being due to 
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a partonic energy loss of the quark that was struck by the virtual photon. 
The colored quark undergoes rescattering inside the nucleus giving rise to 
induced gluon radiation. However, the authors of Refs. 7 and ourselves 8,9 
achieve a very good description of the experimental data by assuming that 
the struck quark forms a color neutral prchadron early after the virtual 
photon-nuclcon interaction. This prehadron then scatters off the surround- 
ing nucleons on its way out of the nucleus. 

In Refs. 8, 9 and references therein we have developed a method to 
incorporate coherence length effects in a semi-classical Boltzmann-Uehling- 
Uhlenbeck (BUU) transport model that uses the Lund string models 
PYTHIA 10 and FRITIOF 11 for the simulation of the elementary DIS pro- 
cess and the (pre)hadronic final state interactions (FSI) respectively. This 
allows for a complete probabilistic coupled-channel description of high- 
energy photo- and electroproduction off complex nuclei. Our simulation 
works on an event-by-event basis and can be directly compared to exper- 
iment accounting for all sorts of kinematic cuts and detector acceptances. 
Our theoretical results are in perfect agreement with the experimental find- 
ings at HERMES. 

2. Model 

In our model we split the lepton-nucleus interaction into two parts: In step 
1) the exchanged virtual photon interacts with a bound nucleon inside the 
nucleus and produces a final state that in step 2) is propagated within 
the transport model. We take into account nuclear effects such as binding 
energies, Fermi motion and Pauli blocking. We also account for shadowing 
of the resolved photon components using the method developed in Refs. 12. 

The final state of the initial photon-nucleon reaction is determined by 
PYTHIA, i.e. the direct and resolved photon-nucleon interactions lead to 
the excitation of one or more hadronic strings which fragment according 
to the Lund fragmentation scheme 13 . Each time a quark-antiquark pair is 
produced in the fragmentation process a string splits into two color neutral 
fragments. If such a color neutral fragment has the mass of a hadron we call 
it a prehadron. After the formation time the hadronic wave-function has 
build up and the prehadron has turned into a hadron with the same energy 
and momentum. In Ref. 14 we have developed a method to extract the 
complete space-time information of each string fragmentation in PYTHIA. 
Although our results should not be overstressed since one applies a semi- 
classical picture to a quantum mechanical problem, our approach makes it 
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Figure 1. Production (top) and formation (bottom) points for a typical Hermes event: 
(photon energy v = 14 GeV, virtuality Q 2 = 2.5 GeV 2 ). The target nucleon is located 
at the origin, the virtual photon is coming in from the left. 

possible to assign the four-dimensional production and formation point to 
each single hadron in each single scattering event. Figure 1 shows the spatial 
distribution of the prehadron production and hadron formation points for 
a typical HERMES event. Obviously, a large fraction of these points fall 
within a distance behind the struck nucleon that is comparable to nuclear 
radii. 

Starting from their production point the prehadrons propagate through 
the nucleus and interact with the surrounding nuclear matter. In this 
work we set their interaction cross section to the corresponding hadronic 
value. The propagation is described within our semi-classical BUU trans- 
port model 9,15 which allows for a probabilistic coupled-channel description 
of the FSI. Each time a particle interacts with a bound nucleon it might 
either scatter elastically or produce new particles that are also propagated. 
At high energies the final state of an inelastic collision is determined by the 
event generator FRITIOF which is also based on the Lund fragmentation 
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scheme. Finally, we end up with a complete lepton-nucleus event that can 
be corrected for experimental cuts and detector acceptance. 



3. Results 

The observable of interest is the so-called multiplicity ratio 

N h {z h ,v) I 
N e (v) \D 

where Nh is the yield of semi- inclusive hadrons in a given (z/ l ,j/)-bin and 
N e the yield of inclusive deep inelastic scattering leptons in the same v-bin. 
The quantity Zh — E^jv denotes the energy fraction of the hadron. For the 
deuterium target, i.e. the nominator of Eq. (1), we simply use the isospin 
averaged results of a proton and a neutron target. Thus in the case of 
deuterium we neglect the FSI of the produced hadrons and also the effect 
of shadowing and Fermi motion. 

Figure 2 shows the charged hadron multiplicity ratio for a 27.6 GeV 
positron beam incident on a 14 N and 84 Kr target. In our calculations 
we assume the same kinematic cuts on the scattered lepton and produced 
hadrons as in the HERMES experiment. In addition we take into account 
the geometrical acceptance of the HERMES detector. The solid curve shows 
the result using the default PYTHIA parameter set. In the calculation 
represented by the dashed curve we used the PYTHIA parameters that 
have been fitted to hydrogen data by the HERMES collaboration 16 . Both 
calculations yield comparable results that are in good agreement with the 
experimental data. 



4. Conclusions 

Our results demonstrate that one needs large prehadronic cross sections to 
describe the experimentally observed hadron attenuation in nuclear DIS at 
HERMES when using the prehadron production points from PYTHIA. In 
this work we have simply set the prehadronic cross sections to the corre- 
sponding hadronic values. So far we have neglected all FSI of the string 
and string fragments prior to the production of the prehadrons and as- 
sumed that string fragmentation in nuclei does not differ from that in vac- 
uum. Both effects might increase the hadron attenuation 7 . In future work 
one might also think of incorporating induced gluon radiation prior to the 
production of the color neutral prehadrons. However, one should prevent 
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Figure 2. Calculated multiplicity ratio of charged hadrons for 14 N and 84 Kr nuclei 
using the four-dimensional production points from the PYTHIA model and setting the 
prehadronic cross section to the hadronic value. The solid line shows the result using 
the default PYTHIA parameters. For the calculation represented by the dashed line the 
fitted parameter set of Ref. 16 has been used. The data are taken from Ref. 3. 

double counting when considering both gluon bremsstrahlung and the FSI 
of the hadronic string. 
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